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ABSTRACT 

The current state of the art in the modehng of pulsar magnetospheres invokes either the vacuum 
or force-free hmits for the magnetospheric plasma. Neither of these limits can simultaneously ac- 
count for both the plasma currents and the accelerating electric fields that are needed to explain the 
morphology and spectra of high-energy emission from pulsars. To better understand the structure of 
such magnetospheres, we combine accelerating fields and force-free solutions by considering models of 
magnetospheres filled with resistive plasma. We formulate Ohm's Law in the minimal velocity fluid 
frame and construct a family of resistive solutions that smoothly bridges the gap between the vacuum 
and the force-free magnetosphere solutions. The spin-down luminosity, open field line potential drop, 
and the fraction of open field lines all transition between the vacuum and force-free values as the 
plasma conductivity varies from zero to infinity. For fixed inclination angle, we find that the spin- 
down luminosity depends linearly on the open field line potential drop. We consider the implications 
of our resistive solutions for the spin down of intermittent pulsars and sub-pulse drift phenomena in 
radio pulsars. 
Subject headings: MHD — pulsars: general — gamma-rays: stars — stars: magnetic fields 



1. INTRODUCTION 

Pulsar magnetospheres are filled with plasma, and the 
presence of plasma affects all magnetospheric properties. 
Although this fact has been well appreciated since the 
early days of pulsar research, the ability to quantitatively 
model plasma effects has emerged mainly in the last 
decade. Generally speaking, the models of pulsar mag- 
netospheres can be classified according to the amount of 
plasma supply they assume. At one extreme is the vac- 
uum magnetosphere, which is modeled with the magnetic 
field of a spinning inclined dipole in vacuum (Deutsch 
19551. As this solution has no plasma, it has no pos- 
sibility of producing any pulsar-like emission. However, 
the fact that this field is described by an analytic for- 
mula has made it the most widely used framework for 
calculating pulsar properties. For example, the charac- 
teristic spin-down energy loss which yielded the fiducial 
pulsar field strength of lO^^G or the polarization sweep 
of th e rotating- vector model ( Radhakrishnan & Cooke 
1969 1 are guided by the vacuum field shape. Slot-gap and 



outer gap models of gam ma-ray emission from pulsars are 
also b ased on this field ( Cheng et al.| 1986a Dyks et al. 
20041. The next order of approximation is models that 
allow plasma emission from the surface of the sta r. These 
include the origina l charge-separated model of Goldre 
ich & Julian (1969 hereafter GJ), and the space charge 



limited mode 



s w ith pair formation (e.g. 



Sutherland 1975 Arons & Scharlemann 



Rud erman & 
9791). i'hese 



models envision both the regions where plasma shorts 
out the accelerating electric fields and the vacuum-like 
regions where acceleration is present (so-called "gaps"). 
This approach allows for more realism in studying plasma 
creation and acceleration, but at the price of being de- 
coupled from the global structure of the magnetosphere, 
as the modification of the vacuum field due to plasma 
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currents is typically not included. Finally, the models 
that concentrate on the global magnetospheric proper- 
ties assume that abundant plasma exists throughout the 
magnetosphere and in the wind. This plasma shorts out 
the accelerating electric fields and provides the corotation 
of field lines with the star. Such models include the rela- 
tivistic magnetohydrodynamic (MHD) description of the 
magnetosphere and its limit for low-inertia magnetically- 
dominated plasmas, the "force-free" description. Which 
of these different regimes is applicable to real pulsars may 
ultimately depend on the physics of plasma supply in the 
magnetosphere . 

Currently, quantitative solutions of the global mag- 
netospheric structure exist only for the vacuum limit 
(Deutsch 1955) and for the lim it of abundant plasma in 
force- fr ee electr odynamics (see Contopoulos et al. 1999 
Gruzinov 2005, |Ti'mokhin|2006| |McKin neyli^006; Komi| 
sarov 200(r] fo r aligned rotators, and Spitkovsky 2006^ 
hereafter S06; |Contopoulos fc Kalap otharakos| |2009| foi 
pulsars with arbitrary inclinations). The real pulsar 
magnetosphere is likely operating somewhere in between 
these limits, with various accelerating gaps, regions of 
pair production, and strong current sheets likely causing 
local violations of the ideal MHD constraint, E ■ B = 0. 
Knowing the structure of the magnetosphere, including 
such non-ideal effects, would be very useful for calculat- 
ing the properties of pulsar emission. Indeed, currently 
the ideal force-free models that include the back-reaction 
of plasma currents on the field structure lack any accel- 
erating fields by construction, and thus cannot be used 
to directly predict the spectra of gamma-ray radiation 
observed by Fermi GST. 

One way to reintroduce accelerating electric fields in 



if 



^ We note that |Komissarov| ( |2006l 
re of the aliened rotator m the reL 



ture of the aligned 
particle inertia. 



has also calculated the struc- 
ativistic MHD limit including 



the magnetosphere is to allow finite resistivity of the 
plasma. Several formulations of resistive force- free equa- 
tions have been pro posed, most notably by |Lyutikov] 
(|2003|) and|GruzinovJ|2007, 2008). In resistive MHD, the 



Uhm s law can be unambiguously defined in the proper 
frame of the fluid by relating the current in that frame 

to the electric field through .7 fluid = o'Enuid , where a is 
plasma conductiv ity (see, e.g. 



Palen- 
owever, the 



Lichnerowicz|1967 
zuela et aL]|2009 ) . In the force-tree system 
fluid velocity along the magnetic field is unknown, and 
only the transverse velocity can be obtained from the 
electromagnetic fields. This introduces some freedom in 
prescribing the Ohm's law in the force - free p icture. In 
the prescription proposed by Lyutikov (2003), the par- 
allel velocity along the field was taken to be zer o. In 

"Strong-Field Electrodynamics" (|Gruzinov| |2007l [2008 



hereafter SFE), the Ohm's law was tormuiated in the 
frame that moves along the field lines with such speed 
that the charge density in that frame vanishes. Despite 
the fact that such a frame formally exists only in space- 
like regions, SFE prescription appears to give smo oth nu- 



merical so lutions throughout the magnetosphere ( Gruzi 
nov 2011), and, most importantly, does not explode in 
current sheets where magnetic fields can go through zero. 
Both formulations are arbitrary, however, because the 
real fluid velocity does not have to follow either frame 
assumption. As an additional constraint, it is useful, 
therefore, to construct a resistive formulation that repro- 
duces physical solutions expected at the extremes of very 
large and very small conductivity of the plasma, namely 
the force-free limit {a -^ oo) and the limit of vacuum 
electromagnetism (a — >■ 0). Below we describe a resis- 
tive prescription that generalizes these schemes and com- 
bines the correct limiting behavior of Lyutikov's scheme 
with current-sheet stability of Gruzinov's formulation. 
We then apply this prescription to numerically calculate 
the structure of resistive magnetospheres in pulsars. 

In the resistive force-free picture of the pulsar mag- 
netosphere, the magnetized neutron star is thought to 
be surrounded by an abundant massless plasma with fi- 
nite conductivity, so that not all accelerating fields are 
shorted out. For simplicity and as a proof of principle, 
we only consider the unrealistic case of constant con- 
ductivity throughout space. More complicated prescrip- 
tions will be studied elsewhere. Our finding is that using 
our formulation of the resistive force-free electrodynam- 
ics we can construct a family of magnetospheres that 
smoothly transition from the Deutsch vacuum solution 
to the ideal force-free magnetosphere as the conductivity 
of the plasma is increased. Such intermediate magne- 
tospheres possess interesting properties that we discuss 
in this paper. We study the variation of the spin-down 
power with magnetic inclination angle as a function of 
plasma conductivity and relate it to physical conditions 
such as the effective potential drop on the open field lines. 
In ^ we discuss the derivation of the resistive force- free 
prescription, in fS] we describe our code for solving the 
equations and present sample magnetospheric solutions. 
Discussion and potential applications to pulsar physics 
are in fJ4] 

2. RESISTIVE ELECTRODYNAMICS 

We describe here our prescription for resistive current 
as used in our numerical code. Ohm's Law is defined 



in the fluid rest frame. In this frame the electric and 
magnetic fields are parallel; otherwise, there would be 
a particle drift across magnetic field lines. The labo- 
ratory frame can be connected to the fluid rest frame 
through two boosts. One boost is in the E x B direc- 
tion in the lab frame and brings the electric field paral- 
lel to the magnetic field. The other boost is along the 
parallel electric and magnetic fields and transforms the 
4-current to the fluid rest frame while leaving the elec- 
tric and magnetic fields parallel. We choose a simple 
Ohm's Law to relate the current and electric field in the 
fluid frame: jfluid = cr-Efluid, where a is the plasma con- 
ductivity. More general formulations of the Ohm's Law 
that take into account time-dependent currents, inertial 
effects, prcssm-e, and the Hall effect have been derived 
(see e.g., Meier 2004). It is likely safe to ignore most of 
these effects m the bulk of the strongly magnetized cold 
pair plasma, as envisioned for typical pulsar magneto- 
spheres; however, singular current sheets may require a 
more elaborate treatment. In this work we use a con- 
stant uniform conductivity throughout the domain and 
restrict ourselves to the simple Ohm's Law to elucidate 
the basic physics. Boosting back to the laboratory frame, 
the current vector can be expressed as 
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See Appendix A for a full derivation of this expression 
and the subsequent limits that we discuss below. Here, 
/3| I is the magnitude of the boost along the parallel elec- 
tric and magnetic flelds to the fluid rest frame, p is the 
charge density in the laboratory frame, and Eq and Bq 
are the magnitudes of the electric and magnetic fields 
in all frames in which they are parallel. Eq and Bq are 
defined by the expressions ( Gruzinov||2007| 2008) 
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where we allow Bq to be positive or negative, depending 
on whether the magnetic field is aligned or antialigned 
with the electric field. ^Not^ that in addition to the advec- 
tion of charge in the E x B direction, there is conduction 
current along both the laboratory frame magnetic field 
and the electric field (the last term in equation flj. There 
is still a principal ambiguity here in defining thenuid rest 
frame, as we have not yet specified the magnitude of the 
parallel boost /3||. In fact, this speed cannot be obtained 
from a purely electrodynamic standpoint without the in- 
clusion of gas dynamics, so a suitable choice has to be 
made. 

One particularly interesting choice of fiuid frame is the 
slowest moving frame that has electric fields parallel to 
magnetic fields, i.e., the frame with /3|| — >■ 0. If parti- 
cles start out on the stellar surface with non-relativistic 
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velocity, then the fluid essentially satisfi es this condition 
in the high conductivity limit (see e.g., Conto poulos & 
Kazanasl 20021 iTchekhovskoy et al.||2009| |tieskini|201 



In this minimal velocity limit the current 



becomes 
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This form of the current has the especially useful prop- 
erty that in the limit of vanishing conductivity the con- 
duction current along the electric and magnetic fields 
goes to zero. Such a magnctosphcre does not develop 
space charge due to the fiow of charge from the star 

along the field lines. Thus, the current in the E x B 
direction also vanishes. This is exactly what we expect 
in the vacuum limit. Since we are interested in producing 
a resistive transition from ideal force-free to vacuum, we 
use t his form of the current in our numerical investiga- 
tion. Lyutikov (2003!) derived an expression for resistive 
current m magnetically dominated plasma, choosing the 
same minimal velocity frame for the formulation of the 
Ohm's law as in our derivation. However, his boost in 
the E X B direction to the fluid frame contains an error 
and does not bring the electric and magnetic fields par- 
allel to one another. As a result, our final expressions 
diS'er. 

We can alternatively express the minimal velocity cur- 
rent as 

j ^ pv + a-Efiuid, (4) 

where the fluid velocity v — c{E x B)/{B^ + Eq) is the 
generalized drift velocity, i^fluid — j{E + vx B), and 7 — 
(1 — v"^ /c^)~^/'^. The drift velocity contains the term E'^ 
in the denominator to account for the non-zero electric 



field in the fiuid frame (see equation A14|. The presence 



of this term allows the current to remain nonsingular in 
current sheets, where the magnetic field vanishes. 
It is instructive to consider other limits for the current. 



Gruzinov ( 2007 2008 1 derived an alternate limit of the 
general current m equation (flj, known as SFE. He postu- 
lated that Ohm's law shoulabe formulated in the frame 
where charge density vanishes, which gives parallel boost 
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where 72 = {B'^ + EI)/{Bq + Eq). We see no special 
reason as to why the fiuid frame charge density must 
vanish. Further, this current does not reduce to vacuum 
as the conductivity drops to zero, as it does in equation 
([3|. Consequently, even for vanishing conductivity the 
SFE solutions resemble ideal non-resistive solutions. 

In the limit of infinitely conductive strongly mag- 
netized plasma, which we refer to as the ideal force- 
free limit, the plasma satisfies transverse force balance, 
pE + j X B ^ 0, and the parallel electric field is shorted 
out, Eq — > 0. However, as cr —)■ co, the product of aEo 
stays finite, and the minimal velocity current becomes 



(Osherovich & Gliner 1988 Gruzinov 1999 Blandford 
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3. NUMERICAL SIMULATIONS 

In order to test the infiuence of global resistivity on the 
structure of pulsar magnetospheres, we employ a three- 
dimensional numerical code (see S06) that imp lements 
the finite differen ce time-domain scheme (FDTD, Taflove 
& Hagness 2005) to evolve electromagnetic flelds from 
Maxwell's equations, 

dtE ^cV X B~ 47rJ, 

dtB = -cV X E, (8) 

where the current is given by equations (l2| and (Is]) . Elec- 
tric and magnetic fields are decentered on the Cartesian 
Yee mesh, as in standard FDTD. Instead of the leapfrog 
time integration typically used in FDTD, we employ a 
third order Runge-Kutta integrator as in S06. This is 
done to remove the interpolation in time that would oth- 
erwise be needed to bring the E and B fields to the same 
time step when computing the source current. We run 
our simulations on a grid of size 1024'^ at Courant num- 
ber 0.5. The central region of our grid is occupied by a 
conducting spherical star of radius R.^, , rotating at angu- 
lar velocity fi, with embedded dipole field of magnetic 
moment jl inclined relative to the rotation axis by an- 
gle a (see S06 for more detail). The electric field inside 
the star is set to ensure rigid corotation, and we smooth 
fields across the stellar boundary as in S06 in order to 
minimize stair-stepping of a sphere on a Cartesian grid. 
The spatial resolution of our simulations is such that the 
light cylinder (LC) radius, i?LC = c/fi, is resolved with 
80 ceUs, and the simulation box is 12.8 i?LC on a side. 
The outer boundary condition is periodic, which limits 
our run time to the light travel time across the grid. We 
set i?* = 30 cells, allowing us to resolve the star while 
keeping the star sufficiently small compared to the light 
cylinder. We run our code for a range of different val- 
ues of the dimensionless conductivity, cr /Q., to produce 
a transition from vacuum at low conductivity to ideal 
force-free solution at high conductivity. We have verified 
that our solutions are converged with spatial resolution, 
as well as run sufficiently long so as to reach a steady 
state in the frame corotating with the pulsar. 

In addition to the resistive solutions obtained with the 
above method, we also consider two special cases of the 
current prescription that we call "force-free" and "vac- 
uum." By force-free solutions we mean "ideal force-free" 
magnetospheres; however, instead of using the full equa- 
tion ([7| for the current, our numerical implementation 
follows the approach in S06, where at each time step 
we evolve the perpendicular part of the current from 
equation (l7| and then correct the electric field by remov- 
ing any electric field component parallel to the magnetic 
field. Further field limiters are used to enforce E < B 
in current sheets. As a result, the force-free scheme of 
S06 has an effective small resistivity, because the cleaning 
of parallel fields can be interpreted as the effect of very 
large conductivity along field lines. Also, as all parallel 



fields are cleaned on every step, such a scheme does not 
have the usual convergence behavior with timestep. The 
resistive scheme using equation (pi gives us a more well- 
defined method for handling magnetospheric dissipation. 
In addition, this resistive scheme does not require special 
limiters to handle current sheetq^ The "vacuum" solu- 
tions described below evolve Maxwell's equations with- 
out volume charges or currents outside the star, allowing 
for significant potential drops along field lines. 

In Fig. Ill we show magnetic field lines in the jl-fl 
plane for 60° inclined force- free dipolc (panel a), resis- 
tive dipoles at (cr/f^)^ = 40 (panel b), (cr/f^)^ = 4 (panel 
c), ia/n)"^ = 0.4 (panel d), {a/n)"^ = 0.04 (panel e), and 
vacuum dipolc (panel f). The field lines are shown after 
1.5 turns, sufficiently long to reach a steady state out to 
several light cylinder radii in the corotating frame. We 
show only the central five light cylinders of our simu- 
lation. Color represents the out-of-plane magnetic field 
into (red) and out of (blue) the page. The maximum 
magnitude of the out-of-plane magnetic field value is dif- 
ferent in each panel, with maximum value increasing with 
conductivity. Table [T] gives the maximum values normal- 
ized to the maximum vacuum valuq3 The maximum val- 
ues occur at the stellar surface and fall on the blue end 
of the color table. To improve contrast, the color table 
shows values up to 30% of the maximum, and a square 
root stretching is applied to the data. Gradients in color 
refiect the sum of in-plane components of conduction and 
displacement currents. 

Fig. la shows the force-free magnetosphere. The gross 
features of this solution are the same as were discussed 



in SO(rl Contopoul os fc Kalapotharak os (2009| and Bai 
[fc Spi tlovsky (20101). I'he polar cap consists of the foot- 
pomts of open field lines that extend out to infinity. Field 
lines with both ends attached to the star form the closed 
field line region, which extends out to a cylindrical dis- 
tance of IRlc- a large scale conduction current flows 
outwards from the polar cap along open field lines and 
returns through the current sheet, current layer, and a 
fraction of the neighboring open field lines. The current 
sheets appear as the sharp transition in the out-of-plane 
magnetic field, starting at the Y-point at the light cylin- 
der. The current layer fiows along the boundary of the 
closed field line region. The conduction current is closed 
in a circuit by a surface current flowing across the po- 
lar cap, a point we will return to when discussing pulsar 
spin-down. 

At (cr/fi)^ = 40, shown in Fig. lb, the qualitative pic- 
ture of the magnetosphere looks quite similar to force- 
free. We find that values ((t/51)^ > 40 give fairly 
good estimates of highly conducting magnetospheres for 
R*/Rlc — 3/8. The dimensionless quantity cr/V, — 
(c/J7)/(c/cr) is the relativistic analog of the Elsasser num- 
ber, defined as the ratio of the Alfven radius in a rotat- 

^ This property is not unique to our resistive current prescrip- 
tion, and is also present in the original formulation by |Gruzinov| 
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20081 



tie discrepancy between maximum field for vacuum and a = 
solutions is numerical and is due to the boundary conditions on the 
star. This is further discussed in Appendix B 

** In comparing to previous force-free results, we note that S06 
rescaled the total magnetic field when displaying the out-of-plane 
component, whereas we rescale only the strength of the out-of- 
plane component. 



TABLE 1 

Maximum out-of-plane magnetic 
field values, normalized to the 

maximum vacuum VALUE, FOR 
EACH PANEL FROM FIG [T] 
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ing system va/^, to the resistive diffusion scale (El- 
1946). The above threshold corresponds to re- 



sasser 



sistive diffusion length scale c/a small compared to the 
light cylinder, c/a < 0.16i?LC ^ ^lCj while satisfying 
c/a < Rpc = i?*(i?*/i?Lc)^^^- We have also run simu- 
lations with our current prescription at higher conduc- 
tivity, {a/fl)'^ = 400, and find spin-down luminosities 
that match the force-free values quite well. Our explicit 
scheme requires very small time step at such high con- 
ductivities though, and i t is perhaps more appro priate to 
use an implicit-explicit (Palenzuela et al. 2009) or fully 
implicit scheme to explore this regime. We note that 
the kink in the current sheet seen in the force-free solu- 
tion near 2i?LC disappears for resistive solutions even at 
high conductivities. It is likely then that this kink is a 
numerical artifact of the cleaning procedure used in our 
force-free formulation. 

Displacement currents play an important role in the 
highly conducting magnetospheres of Figs, la and lb. 
The displacement currents in the current sheets actu- 
ally dominate over conduction currents for the 60° rota- 
tor. Displacement currents are not present in the mag- 
netospheres of aligned rotators, but the strength of the 
displacement current in the current sheets of highly con- 
ducting magnetospheres increases monotonically with in- 
creasing inclination angle. This is in contrast to the con- 
duction currents, whose strength is independent of incli- 
nation angle for fixed conductivity. Both displacement 
and conduction currents weaken considerably with de- 
creasing conductivity, reflected by smaller gradients in 
the out-of-plane magnetic field in panels (b) through (e) 
of Fig.jl] Recall that the maximum in the color table cor- 
responds to weaker out-of-plane magnetic field at lower 
conductivities. The extent of the closed field line re- 
gion also expands with decreasing conductivity, evident 
in the same sequence of panels. We further expect the 
current sheet to thicken with decreasing conductivity. At 
(a/n)'^ = 0.04, shown in Fig. le, the current sheet should 
have characteristic width c/a of order BRlc- 

There is an obvious increase in the radius of toroidal 
field sign reversal beyond the light cylinder as one tran- 
sitions from the force-free to the vacuum magnetosphere. 
This can be understood from the increase in the winding 
radius of the spiral pattern for lower conductivities. The 
winding of the spiral is determined by the characteristic 
speed of the radial outfiow multiplied by the pulsar pe- 
riod. In the force- free case, the radial outflow is given 
by the radial component of the E x B velocity, which in- 
creases with radius so that the Lorentz factor grows lin- 
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Fig. 1. — Magnetic field lines in the fl — Q plane for a 60° inclined dipole. Color represents out-of-plane magnetic field into (red) and 
out of (blue) the page. The color table shows only values up to 30% of the maximum of the out-of-plane magnetic field, and a square 
root stretching has been applied to its magnitude. The maximum out-of-plane magnetic field values are given in Table IT] Conduction and 
displacement currents weaken with decreasing conductivity, (a) force-free dipole; (b) resistive dipole at (a/fl)'^ = 40; (c) (a/fl)'^ = 4; (d) 
(a/n)'^ = 0.4; (e) {cr/nf = 0.04; (f) vacuum dipole. 
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early with cylindrical radius: jexb = a/I + (R/Rlc) 
{ Contopoulos & Kazanas 2002 1. Thus, the force- free out- 
tiow approaches the speed ot light only asymptotically. 
In the vacuum case, the "outflow" is always at the speed 
of light, so we expect the sign transition of the field to 
occur near the half-wavelength of the vacuum wave, or 
7ri?LC, as is confirmed by Fig. If. The force-free case 
of Fig. la reverses near 2i?LCj indicative of the smaller 
outflow speed near the light cylinder. 

The surface currents flowing across the polar cap of 
the pulsar exert a spin-down torque on the star. The 
rotational energy loss is connected to infinity by an out- 
ward directed Poynting flux, which can be thought of as 
the ejection of toroidal field. We compute the Poynt- 
ing flux at the light cylinder for each of our simulation 
runs and show the results in Fig. [2] All spin-down val- 
ues are calculated for i?» = 3/8i?LC- The spin-down 
curves have been normalized to Lq, defined as 3/2 times 
the power of the orthogonal vacuum rotator with fi- 
nite i?* = 3/8Rlc- Two-dimensional axisymmetric cal- 
culations for small star tell us that the spin-down of 
the aligned force-free rotator should be 3/2 times larger 
than the power of the orthogonal point vacuum dipole. 



L^_=2j?£i^/3^ 



Contopoulos et al. 1999 



^ ^Gruzinov] 

2005||McKinney|20(]6| S06J. I'he short dashed line shows 
the function L/Lq = 1 + sin a, the expected force-free 
curve in the limit of small star with inclined dipole (S06) . 
The gray band around the force-free curve indicates the 
uncertainty in the measurement due to boundary effects 
and numerical dissipation of Poynting fiux in the magne- 
tosphere (see Appendix B for thorough discussion) . Sim- 
ilar bands are implied but not shown for resistive spin- 
down curves. The long dashe d line shows the an alytic 
vacuum Deutsch field solution ( Michel fc Li| |1999). One 
of our principal results is that we see a smooth mono- 
tonic transition from the force-free spin-down curve to 
the vacuum spin-down curve for decreasing conductivity. 
All spin-down curves in Fig. [2] show a strong depen- 
dence on both inclination angle and conductivity. The 
radial Poynting flux carrying the spin-down power is pro- 
portional to the product of poloidal and toroidal mag- 
netic fields {EqB^ ~ URBpB^/c), both of which are af- 
fected by the strength of displacement and conduction 
currents in the magnetosphere. One can get a qualita- 
tive feel for the relative significance of two contributions 
by considering the limiting cases. In vacuum, the in- 
crease of spin-down with inclination is solely due to ris- 
ing displacement current. This current is likely respon- 
sible for much of the angular dependence of the resistive 
and force-free solutions. The increase in spin-down with 
increasing conductivity is due to the additional sweep- 
back and opening of the poloidal field brought on by the 
increasing conduction current. We find that spin-down 
luminosity at any conductivity has an angle-dependent 
term proportional to sin a. We parameterize the spin- 
down curves in Fig. ^with functions /(cr/fi) and g{a/Vl) 
such that L/Lq = f{cr/^) + g{<j/^) s\v? a. We find the 
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Fig. 2. — Spin-down luminosity dependence on inclination angle 
for force-free, a sequence of resistive, and vacuum dipoles. Spin- 
down is normalized by 3/2 times the spin-down power of the or- 
thogonal vacuum rotator. We see a smooth monotonic transition 
from force-free to vacuum with decreasing conductivity. 



piecewise linear fit 
L 

L 

0.004 < (a/n)^ < 0.4 



= 0.3 + 0.3 log{a/nf + 1.2 sin^ a , (a/n)^ > 0.4; 

= 0.2 + 0.08 log(cr/f7)2 + (1.3 + 0.2 log{a/n)^) sin^ a, 

(9) 



The amplitude of the angular dependence, g{a/Q), is 
constant for (cr/J7)^ > 0.4 and begins to transition to 
the vacuum value below (a/H.)'^ = 0.4. 

We have thus far shown how spin-down luminosity de- 
pends on plasma conductivity, but the physical meaning 
of the conductivity is not entirely clear. It is instructive 
to reinterpret the conductivity parameter, a/fl, in terms 
of the potential drop along open field lines in the coro- 
tating frame. This gives us a handle on the deviation 
of the magnetosphere from ideal force-free, which has 
vanishing potential drops along field lines. The electro- 
magnetic fields in the frame corotating with the pulsar 
are obtained via a coordinate transformation from the 
laboratory frame ( |Schiff|p39l [Gr0n 1984): 



and 



E' = E+ X B 



B' ^ B. 



(10) 
(11) 



Since the fields are steady in the corotating frame, V x 
E' = Q and the corotating electric field can be written as 
the gradient of a scalar potential, i.e., E' — Vx- Taking 
the line integral of the corotating electric field along a 
magnetic field line, Z, we find 



Ax 



E' -dl 



E-dl = Vdr. 



(12) 
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We see that the potential drop along field lines in the 
corotating frame can be computed directly from the lab- 
oratory frame fields. Although in resistive solutions par- 
ticles will actually drift across magnetic field lines, in ad- 
dition to accelerating along them, we choose to study the 
field-aligned potential drop Vdrop SiS a fiducial measure of 
particle energy gain. In fact, the exact particle trajec- 
tories we choose make little difference when estimating 
potential drops because the electric field is potential in 
the frame corotating with the pulsar. 

Consider field lines starting on the stellar surface in 

the /2— $7 plane separated by 15° in latitude from pole to 
equator. For every such field line we integrate the electric 
field to find the maximum potential drop along each field 
line. We then determine the field line with the largest 
overall potential drop for a given magnetic inclination. 
Integrating field lines separated by 15° increments on 
the stellar surface is sufficient to give a good estimate 
of the maximum potential drop along field lines. Fig. 3] 
shows the maximal potential drop as a function of dipole 
inclination angle for different conductivities. All results 
have been normalized to the potential drop from the pole 
to the equator of an aligned vacuum rotator in the lab- 
oratory frame, Vb — \fl\/{RhcR*)- As our models do 
not prescribe a high conductivity to the closed field line 
region, the available accelerating potential is generally 
limited by the pole-to-equator potential difference, rather 
than the smaller polar cap potential, Vpc = VoR^/Rlc, 
with a few notable exceptions. Low conductivity solu- 
tions, (cr/r2)^ < 0.04, at high inclination angle, a > 45°, 
have Vdrop scaling intermediate between Vb and Vpc . The 
orthogonal rotator drop at (cr/fi)^ = 0.04 scales roughly 
with Vpc. Using the vacuum Deutsch fields, one can 
obtain that the aligned vacuum rotator potential drop 
scales exactly proportional to Vb, whereas the orthogonal 
vacuum rotator drop scales closer to Vpc. From Fig. [3] we 
see that for (cr/ri)^ > 0.04 the potential drop is roughly 
independent of inclination angle for fixed cr, i.e., there 
is a one-to-one map between conductivity and maximal 
potential drop along a field line originating in the fl ~ Q 
plane for all inclination angles. The maximal potential 
drop increases from zero for force-free (or a/ ft — > cx)) to 
Vdrop/Vo ~ 0.2 for (a/n)^ = 0.04. 

The fact that potential drop is independent of inclina- 
tion angle for (u/ri)^ > 0.04 allows us to relate luminos- 
ity to potential drop in a very simple manner. We present 
the result here and explain its origin in more detail be- 
low. Fig. [4] shows the spin-down luminosity as a function 
of potential drop for inclination angles a = 0, 30, 60, 90°. 
Spin-down luminosity increases with increasing inclina- 
tion angle and with decreasing field line potential drops. 
We fit the spin-down curves in Fig. |4] with the linear 
relation 



L 



0.9 1 



Vdr. 



0.2Vb 



l.lsin^a, 0< Vdrop <0.2Vo- 

. (13) 

This formula applies only to the domain for which we 
have data shown with solid lines in Fig. [4| i.e., < 
Vdrop /Vb < 0.2. A similar spin-down fo r mula was con- 
structed by Contopoulos fc Spitkovsky] ( |2006| for the 
case of a finite gap at the base of the open field lines 
in a force-free magnetosphere. The spin-down was also 




((7/n)'=40 



20 40 60 

Inclination angle, degrees 



Fig. 3. — Open field line potential drop with inclination angle for 
force-free, a sequence of resistive, and vacuum dipoles. Results are 
normalized to Vq = \fi\/Ri^cR*- Fo'' (o"/^)^ > 0-04 the potential 
drop is relatively independent of inclination angle. 
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Fig. 4. — Spin down dependence on the open field line potential 
drop for a sequence of different inclination angles, a = 0° , 30° , 60° , 
and 90°. Potential drop is normalized to Vb = Ifil/R'LcR* a^nd 
luminosity to 3/2 the spin-down of the orthogonal vacuum rotator. 
For Vdrop/M) < 0.2 spin-down falls linearly with potential drop. 

found to be linear in the potential drop in the gap but 
had a different slope and angular dependence of the form 
(1 — Vdrop/Vpc) cos^ a + sin^ a. This angular dependence 
was not a rigorous derivation, but a prediction based 
on aligned force-free and orthogonal vacuum limits. We 
note that our spin-down formula reduces to the force- 
free spin-down formula from S06 within error bars when 

^drop/K) = 0. 

For Vdrop/Vb > 0.2 the dependence of spin-down on 
the potential drop enters a different regime (see dotted 



lines in Fig. [4]). We do not provide a fit here, as po- 
tential drop no longer scales with Vq at high inclination 
angle. The potential drop along field lines in this regime 
is due to the reemergence of the vacuum electric fields as 
conductivity is reduced. There are three contributions 
to the vacuum electric field which are provided by the 
quadrupolar surface charge on a finite-radius star and 
by central and surface monopolar charges. The central 
monopolar component provides part of the radial electric 
field needed for corotation of the magnetized stellar in- 
terior. If the star initially is uncharged (as is assumed in 
our sinuilations), a surface charge of opposite sign com- 
pensates the net charge of the interior of the star (see 
Michel fc Li||1999 for thorough discussion). This surface 
charge, together with the induced quadrupolar surface 
charge from corotation, can leave the star and be redis- 
tributed throughout the magnetospherc when conductiv- 
ity is turned on, or when the work function of the surface 
is low. Such redistribution lowers the maximum available 
potential drop. 

The available mo nopolar surface cha rge varies with in- 
clination as cos a ( [Michel fc Li 1999), so the potential 
drop associated with it disappears tor orthogonal rota- 
tors. Hence, the value Vd^op ~ 0.2Vo for the vacuum 90° 
rotator in Fig. [3] represents the quadrupolar contribution 
to the potential drop. At conductivity (a/fl)^ ^ 0.04, 
the maximal potential drop is close to the orthogonal 
vacuum rotator value of Vdrop ~ 0.2Vo at all inclina- 
tion angles. Qualitatively, magnetospheres with finite 
conductivity {a/i^)'^ > 0.04 can effectively redistribute 
the surface charges and maintain the linear relationship 
between the maximal potential drop and the spin-down 
luminosity. Below (a/fi)^ = 0.04, the surface charge re- 
distribution is incomplete (the monopolar surface charge 
is concentrated closer to the star in the corotating steady- 
state solution) and the curves in Fig. p^ gain a tilt more 
reminiscent of the vacuum solution. Tne potential drops 
are not constant with inclination ang le at a given con- 
ductivity, and the linear relation ( 13 1 no longer holds. 
We have explicitly verified that by the end of the simula- 
tion these solutions reached steady-state in the corotat- 
ing frame. 

Redistribution of surface charge in near- vacuum mag- 
netospheres also affects the spin-down power. Look- 
ing back at Fig. [2] and the piecewise linear relations 
(l9f , we see that the final 30% of the spin-down tran- 
sition between force-free and vacuum rotators occurs for 

{a/n)'^ < 0.04. Despite the fact that potential is no pulsars switch between an 

longer constant with angle for each value of the conduc- 
tivity in this regime, the transition in spin-down is still 
continuous and smooth. 

4. DISCUSSION 

We have presented a continuous one-parameter family 
of pulsar magnetospherc solutions that span the range 
between the vacuum and force-free limits. Each solution 
of the family is characterized by the value of the con- 
ductivity parameter, a/Q, which is related to the max- 
imum potential drop that a test particle can experience 
as it moves along field lines, Vdrop- The zero conductiv- 
ity limit, a/fl — > 0, Vdrop — > Vdrop, vac, yiclds the vacuum 
magnetospherc. It shows a substantial effective potential 
drop that makes up a significant fraction of the rotation- 
induced potential difference between the pole and the 



equator of the star (20—60% depending on the inclina- 
tion). The high conductivity limit, a/il — > oo, yields 
essentially an ideal force-free magnetospherc. Unlike in 
vacuum, abundant magnetospheric charges short out any 
potential differences along field lines, leading to vanish- 
ing effective potential drop, Vdrop, ff -^ 0. While in the 
vacuum case nearly all field lines, even those that ex- 
tend beyond the light cylinder, eventually return to the 
star and are formally closed, in the ideal force-free case 
a fraction of field lines open up and reach infinity. Our 
simulations show that resistive high-cr pulsars spin down 
at least 3 times faster than resistive low-cr and vacuum 
pulsars (for the same value of inclination), in agreement 
with earlier ideal force- free simulations (806). Our resis- 
tive solutions bridge the gap between the force-free and 
vacuum limits by having intermediate values of major 
magnetospheric parameters, such as effective potential 
drop, < Vdrop < Vdrop. vac J the fraction of open field 
lines, and the spin-down rate. 

Before discussing possible observational implications of 
these solutions, we note that our approach to modeling 
resistive pulsar magnetospheres is quite simplistic: we 
use a form of Ohm's Law in which we neglect several 
terms ( Meier 2004 ) and we assume constant conductivity 
a/fl throughout all space. The finite bulk conductivity 
used in this work can be thought of as not due to col- 
lisions, but due to the insufficient plasma supply in the 
magnetospherc that cannot short out all the accelerat- 
ing fields. It is possible that real magnetospheres have a 
lower (anomalous) conductivity in the current layer and 
the current sheet than in the rest of the magnetospherc. 
We might also expect a lower value of conductivity on the 
open field lines than on the closed field lines. Since our 
models do not explicitly model the plasma fluid (Q, the 
plasma velocity component along the direction of mag- 
netic field is unconstrained, which necessitates the choice 
of a frame in which to write the Ohm's Law. We choose 
the so-called minimum velocity frame, which has a num- 
ber of attractive properties (see Spl), yet this is not a 
covariant choice. With these important caveats in mind, 
let us now for the sake of argument take our model seri- 
ously as a description of pulsar magnctosphere and con- 
sider the consequences. 

Intermittent pulsars offer a uni que test bed of pul 



2006 



sar theory (e.g., PSR B1931-F24 [Kramer et al.|, 

and PSR J1832-h0029, |Lyne||2009| se e also Wang et al'.l 

2007[ [Zhang et aL[|2007( and ITimokhin 2010). Such 



on" , radio-loud, state in 
which they behave like normal radio pulsars, and an 
"off", radio-quiet, state in which they produce no de- 
tectable radio emission. The two intermittent pulsars, 
for which published data exist, have quite different duty 
cycles: PSR B1931-t-24, with a period P sa 0.8 s, cy- 
cles through the "on" -"off" sequence of states approxi- 
mately once a month, whereas PSR J1832-I-0029, with 
a period P ~ 0.5 s, kept quiet for nearly two years. 
The spin-down rate for each of these pulsars is larger 
in the "on" state than in the "off" state by a factor 
/on-foff — 1.5. Such a substantial difference in spin-down 
rates suggests that the pulsar magnetospherc undergoes 
a dramatic reconfiguration as it transitions between the 
"on" and "off" states, yet such a transition was reported 
for PSR B1931+24 to take place in just over 10 pulsar 
periods. 
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Kramer et al. (2006) propose that in the "on" state, 
plasma tiUs the pulsar magnetosphere and supports 
(poorly understood) plasma processes that produce ra- 
dio emission. In this picture, the pulsar transitions to the 
"off" state when, due to some unknown trigger, magneto- 
spheric pair production shuts off: the remaining plasma 
leaks off the open field lines and the pulsar goes radio- 
quiet. In this scenario, we expect the pulsar in the "on" 
state to spin down due to force-free energy losses and in 
the "off" state due to vacuum dipole losses. However, 
this picture is ruled out by observations: the force- free 
spin-down rate is larger than the vacuum spin-down rate 
by a factor, fg^vac > 3 (S06), that is clearly inc ompat 
ible with the obse rved value, /on^off — 1-5 (e.g., Beskin 
fc Nokhrina||2007l ) 



Uur resistive magnetospheres provide a possible res- 
olution to this problem since they have intermediate 
spin-down rates between vacuum and force-free. Let 
us associate the "on" state with the force-free mag- 
netosphere and the "off" state with one of our inter- 
mediate resistive magnetospheres. According to equa- 
tion ( |13[ ), for the ratio of spin-down powers in these 
two states to equal /on->-off; the potential drop has to 
be Vdrop = M)(/o„^off - 1) (0.9 + l.lsin2a)/4.5/o„_>off. 
For /on->-off ~ 1-5, the ratio Vdrop/Vo increases monoton- 
ically with increasing inclination angle, a, from 0.07, for 
an aligned rotator (a = 0°), to 0.15, for an orthogonal ro- 
tator (a = 90°). Plugging in for Vq, using the measured 
spin-down parameters of PSR B1931-I-24 (we associate 
the observed spin-down rate in the "on" state with the 
spin-down rate of the force-free magnetosphere) , and as- 
suming a neutron star mass, M* = IAAIq, we obtain 
Vdrop = 1-5 X 1016(0.9 + l.lsin^a)!/^ [y]. This value 
significantly exceeds the characteristic val ue of the po- 
lar cap potential drop exp ected for pulsars ( Goldreich & 



Julian||1969 Arons 2009), so the pulsar should have no 
difficulty m creating pairs and shorting out the large po- 
tential. The reason for this discrepancy is our simplifying 
assumption that the conductivity is constant throughout 
the magnetosphere. This assumption causes the poten- 
tial drop not just across the polar cap but across the 
whole stellar surface. In real pulsars, we expect that the 
closed field line region is filled with plasma and is highly 
conductive. This plasma shorts out parallel electric field 
in the closed field line region and screens the potential 
drop everywhere except across the polar cap. The poten- 
tial drop across the polar cap, Vpc, is « {R*/Rlc) times 
the pole-to-equator potential drop, Vq, which gives 
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V, 
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on— ^ off 
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4.5/on-)-off 

3.8 X 10^2 [V], 



(14) 



for PSR B1931+24, where we assumed R^ = 10 km. 

Addition of resistivity can, in principle, allow the de- 
tailed study of pulsar spin-down and the determination 
of pulsar braking indices. For a spin-down law fl = Q{t), 
the braking index is defined as n = flCl/Cl^. In order 
to measure a braking index numerically, we run a series 
of simulations for different values of angular frequency, 
n, and map out the dependence, L{U). We find that 
L{fl) ex fl'^ for both purely force-free and vacuum solu- 
tions, which translates into a value of the braking index. 



n = 3, as we will see below. This i s larger than the 
range of observed values, n ~ 2.5—2.8 (Livingstone et al. 
20071. One way to lower the braking index is to allow 
the evolution of the extent of the magnetosphere (the ra- 
dius of the Y-point) to lag behind the outwar d expansion 



of the light cylinder du e to spin down (e.g., Contopou- 



los & Spitkovsky 2006). The extent of the Y-point is 
controlled by the rate of reconnection at the edge of the 
magnetosphere. We calculated the braking index in our 
simulations by comparing the magnetosphere with con- 
stant conductivity spun up to different periods, and we 
do not find significant deviations from n = 3. This rapid 
spinup biases our answer, however, because the Y-point 
extends to the light cylinder within one rotation period. 
It is not entirely clear that we would obtain the same 
braking index if we allowed the pulsar period to increase 
during a single simulation as the pulsar spins down and 
the light cylinder slowly recedes. The extent of the Y- 
point may depend on the past history of the pulsar. To 
do this calculation, we need to increase the pulsar period 
several times within a single simulation and measure the 
resulting spin-down luminosity after each increase. We 
must also ensure that the reconnection at the Y-point is 
controlled by the resistivity we impose and not by nu- 
merical effects. A more accurate braking index calcula- 
tion would require much longer simulations, and possibly 
higher spatial resolution near the Y-region. 

Another possibility that can affect the braking index 
is the evolution of conductivity with time. Such time- 
evolution translates into evolution in pulsar dimension- 
less luminosity, £ ~ L/Lq (see Fig.[2|, and leads to n =^ 3: 



an 



3-2 



d£/dt 



(15) 



where t = P/2P is the pulsar age, and we used the 
scaling, Lq cx il*. If dimensionless plasma conductiv- 
ity does not change in time, ct/H — constant, then £ = 
constant and equation ( 15 1 gives n = 3. Purely force-free 



and purely vacuum pulsars, discussed above, fall into this 
category. As a pulsar ages, it is natural to expect that its 
plasma supply depletes and the pulsar becomes progres- 
sively more vacuum-like: in other words, both cr/f2 and £ 
decrease in time (see Fig. [2]). For such pulsars, d^/di < 0, 
and equation ( 15 ) gives n > 3. Of course, the converse 
is also true, and the braking index lower than 3 would 
result if a pulsar becomes progressively more force-free- 
like, i.e., if pulsar dimensionless luminosity increases in 
time, d£/At > 0. Presently, however, it is not clear why 
such behavior would be physically expected. Thus, the 
addition of bulk resistivity does not seem to easily solve 
the braking index puzzle, and a more elaborate explana- 
tion is still required. 

Our resistive magnetospheres naturally lend them- 
selves to modeling the sub-pulse drif t phenomena ob- 
served in a number of pulsars (see e.g.,IAskeg ar fc Desh 
pandc 2001; Deshpand e fc Ran kin 2001; Rankin et al 
p03| [Contopoulos fc Spitkovsky] |2006| ). The following' 
discussion takes place in the corotatmg frame, where we 
define electromagnetic fields by equations (10) and (111. 
In this frame there is transverse particle dritt with veloc- 
ity of order Udrift ~ E± x B/B^. This drift velocity can 
be thought of as the minimal velocity of particles, ignor- 
ing motion along the magnetic field lines. The transverse 
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electric field E± reverses sign through the axis of maxi- 
mum potential drop, and so there is rotation of field lines 
about this axis. In general, this axis is not in the direc- 
tion of the magnetic moment, and in fact does not even 
need to be a straight line. If we consider a simplified pic- 
ture of the magnetosphere in which there is no parallel 
electric field in the closed field line region, the potential 
drop across open field lines is given roughly by the longi- 
tudinal potential drop along field lines. The differential 
rotation can then be roughly related to the longitudinal 
potential drop by 

An = AVd,op/{BpRl) = Al/drop/*cap. (16) 

In the laboratory frame, the field lines are undergo- 
ing rotation about the pulsar spin axis (with the angular 
frequency il), and in addition they are rotating differ- 
entially around the maximum potential axis in a retro- 
grade sense. The differential rotation of plasma can ex- 
plain sub-pulse emission features that drift with respect 
to the otherwise periodic light curve of the pulsar. Our 
result is in contrast to the standard pictur e in which the 



plasma precesses about the magnetic axis ( Ruderman & 
Sutherland] |1975|). The magnetic colatitude parameter 



in the cartogra pfiic transformation of, e.g. |Deshpande fc| 
[Rankin' (2001^ , must then be modified accordingly. 

To better visualize the differential rotation in the con- 
text of our resistive magnetospheres, consider the full 
drift velocity 

E X B 

Vdntt = r.9 , ^2 *^' 



B^ + E^ 



(17) 



again in the corotating frame. Fig.lSlshows magnetic field 
lines in the /2 — il plane for a 60° inclined dipole at con- 
ductivity (cr/fi)^ — 4. Color indicates out-of-plane drift 
velocity, with red (blue) representing velocity into (out 
of) the page. The velocity has been rescaled by raising 
its magnitude to the 1/2 power. The axis of maximum 
potential extends from the star roughly halfway between 
the rotation and magnetic axes. The differential rotation 
of plasma about this axis is denoted by the transition 
from red to blue across the axis, indicating a reversal in 
the sign of the out-of-plane drift velocity. Note that the 
axis bends off in the direction of open field lines beyond 
the light cylinder. These corotating frame drift velocities 
can be quite substantial. The maximum illustrated drift 
velocity occurs near the current sheets and corresponds 
to a velocity of 0.8c. The differential rotation is much 
stronger than that implied by sub-pulse drift, but we 
expect the qualitative features of differential rotation to 
persist in solutions with more realistic open field line po- 
tential drops. At lower conductivities the magnitude of 
differential rotation about the axis of maximum potential 
drop increases, and the conductivity is a parameter that 
can be tuned to attempt to match observed sub-pulse 
emission features. 

Incidentally, the bundle of field lines surrounding the 
axis of the maximum potential corresponds to the field 
lines that carry the largest di stributed curre nt density 
on the polar cap (see Fig. 4 in Bai & Spitkovsky 2010). 
If the strength of the current is associated with radio 
emission, then the core emission would have a centroid 
that is offset from the magnetic pole. The offset could 
be as large as half of the polar cap radius. This would 




Fig. 5. — Magnetic field lines in the ft — U plane for a 60° in- 
clined dipole at conductivity {cr/f!)-^ = 4. Color represents the 
out-of-plane drift velocity in the corotating frame, with blue color 
representing the out of the page direction. There is retrograde 
differential rotation about the axis with maximum potential drop, 
situated roughly halfway between the rotation and magnetic axes 
when inside the light cylinder, as well as about the boundary of 
the closed field line region. The maximum value on the color table 
corresponds to a drift velocity of 0.8c. 

introduce potentially significant modifications to the po- 
larization sweep of the core radio emis sion and cause de- 
viations from the expected S-curve of Radhakrishnan & 
Cooke (1969). For a number of pulsars, it would imply 
differences m the inclination and viewing angles inferred 
from the shape of the polarization sweep. 

We conclude with prospects for future research. We 
will generalize our assumption that a/fl is constant in 
the magnetosphere. In addition to having higher con- 
ductivity in the closed zone, we will experiment with pre- 
scriptions for anomalous resistivity in the current sheet. 
These improvements will yield a more realistic magneto- 
spheric structure, that can be accurate enough for geo- 
metrical modeling of gamma-ray light curves from pul- 
sars. The gamma-ray pulse formation is sensitive to the 
geometry of the magnetosphere, and , in particular, to 
the fi eld lines near the current sheet (Bai & Spitkovsky 
20101. Deviations from the force-free geometry may be 



important for modeling the light curves of older p ulsars 
which require w ider gaps in the outer-gap models (Wat- 
ters et al.||2009| 
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APPENDIX 

A. RESISTIVE CURRENT DERIVATION 

We derive here in full generality our resistive current prescription. We start in the fluid frame, with charge po and 
current flowing along the electric field with magnitude aEg. For convenience, we pick the current and electric field 
to point along the positive z axis. The magnetic field can point along the positive or negative z axis, depending on 
whether Bq is positive or negative. If we boost along the z axis with /3i — (0, 0,l3z), the current 4 vector in the boosted 
frame satisfies 

Iz -Pzlz 



pc 



z'~iz 



Iz 



Pac 




aEo 



(Al) 



Ix 
Pxlx 




X 


\p'c^ 







1_ 




[fz\ 



Boosting again in the x direction transverse to the electric and magnetic fields with /32 = {Px, 0, 0), 

pc 

Jx 

iz 

we obtain for lab frame quantities the system of equations 

P = IxlzPo - lxl3zlz<yEo/c, 

jx = -fixcp, 
jz = -Pzlzcpo + ^zO-Eq. 



Rearranging equation ( A3 ) 



Po 



Ixlz 



- (S^aEo/c, 



(A2) 

(A3) 
(A4) 
(A5) 

(A6) 



12 

and plugging into equation ( A5 ) , 

jz = -PzlzCi (- PzO-Eq/c) + "fzdEo 

Ixlz 

= -^cp+—aEQ. (A7) 

^x "tz 

We next determine the magnitude of the perpendicular boost in the x direction. The electromagnetic fields are 
invariant under the boost /3i . After applying the boost /32 , the lab frame electromagnetic fields become 

-B, =0 Ey = -i,[E'y - p,B',) = ~-i.,p.,Bo E, = -f,{E', + (3,B'y) ^ j,Eo 

B^^O By = -i,[B'y + p.,E',) = -f^p.,Eo B, = -f,{B'^ + p.,E'y) = j,Bo. (A8) 



Hence 



B' = ^'f°t5 ^ (A9) 



/3. = 15^^- (AlO) 



and 
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-Bi 
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+Er 


B^ 


+ El 



ll = i^-^- (All) 



B',+Et, 



The lab frame x direction corresponds to the —E x B direction. The y and z coordinates are rotated with respect to 
the B and E±^ directions. We let 6 denote the angle between the z axis and the laboratory frame magnetic field and 
assume without loss of generality that jS^ > 0. The current in the E x B direction is then 



Noting that 




E^oB^o/B^ (A13) 

(A14) 

The boost — /3i = {—/3x,0,0) is the slowest boost from the lab frame that brings the electric and magnetic fields 
parallel to one another, and /Jj. can be thought of as a generalized drift velocity. From (3^^ we immediately obtain the 
transverse current 

The current in the direction of the magnetic field is 

3s ^3z cose ^ ( pc+ —aEo)——- 

B Ix Iz B 

( R ^> ^. BqB + EqE-B/B 
= {-P,pc+—aEo) ^2 ^^2 ■ (A16) 

The current in the direction E±^ is 

■ ■ a I P^ , 1 p \lxl3xEo 
J-^ = JzSva^O = ( pc^ uEq) — 

^i Ix Iz B 

H-&.c+|.E„,Jf±,^ (A17) 

Hence the laboratory frame current vector can be expressed in full generality as 



^ pcExB + {-p,pc + V|^(l - Pl)aEo){BoB + E„E) 

^- ^-^^^ ■ (^1«) 

In ^f3z has been replaced by /3|| to emphasize that it is the magnitude of the boost along the parallel electric and 
magnetic fields to the fluid rest frame. 
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Fig. 6. — Poynting flux with radius for force-free, a sequence of resistive, and vacuum dipoles inclined at o = 60°. The resistive Poynting 
flux values slope dow^nwards with increasing radius, reflecting magnetospheric dissipation. Results are normalized to the force-free spin-down 
value on the star. 

B. NUMERICAL EFFECTS ON SPIN-DOWN POWER 

In principle we should be measuring spin-down values at the stellar surface, but there is uncertainty in the stellar spin- 
down measurement due to stair-stepping at the spherical inner boundary on our Cartesian grid. To avoid this issue, we 
measure spin-down at the light cylinder. Dissipation of Poynting flux inside the light cylinder artificially suppresses our 
spin-down estimates measured at the light cylinder, though. We quantify these uncertainties by computing the force- 
free spin-down values measured both on the star and at the light cylinder. The true force-free spin-down luminosity 
likely falls within the shaded grey region in Fig. [2] bounded by the stellar and light cylinder spin-down values. There 
is also a deviation from the vacuum Deutsch solution for our zero conductivity solution (compare dashed line to 
(a/n)'^ = line). This difference is due to our boundary conditions at the star. The smoothing of the fields across 
the stair-stepped boundary leads to small charge bleed-off from the inside of the star to the outside. The first term in 
equation pi) is then not identically zero just outside the star even in the a = case. This numerical artifact does not 
affect the solutions with high conductivity when the physical current exceeds the numerical smoothing current, but 
below (cr/fi)^ = 4 X 10"'^ it can influence the field structure and slightly modify the spin-down power. 

To better understand these numerical issues it is helpful to look at the run of Poynting flux with radius. In Fig. [6] we 
show the Poynting flux integrated over spherical shells of varying radii for force-free, vacuum and a range of resistive 
solutions with a = 60°. The results are normalized to the force-free value at the star, located at R^, = 0.375Rlc: 
and we show the run of Poynting flux with spherical radius out to 2.25Rlc- The vacuum and (a/ri)^ = curves are 
flat, indicating negligible dissipation of Poynting flux with increasing radius, but there is an offset between the curves 
that we attribute to our imperfect inner boundary. Force-free is in principle dissipationless and should have the run 
of Poynting flux flat with radius. As was argued above, our method for cleaning parallel electric field in force-free 
simulations is resistive and causes Poynting fiux to slope downwards with increasing radius. Inside the light cylinder, 
the drop is due to artificial volume j • E dissipation above the polar caps. This dissipation varies with magnetic 
inclination angle and is responsible for the varying width of the gray error band with angle for force-free solutions in 
Fig. [2] Outside the light cylinder, the drop in Poynting flux is due to it disappearing into the current sheets. 

The resistive runs in Fig. [6] have Poynting flux dropping with increasing radius, as we expect, indicating dissipation 
throughout the magnetosphere. Inside the light cylinder the dissipation is primarily physical j ■ E dissipation above 
the polar caps. There is also some dissipation due to the disappearance of Poynting flux into the current layer, though 
signiflcantly weaker in magnitude as compared to the volume j ■ E dissipation. Outside the light cylinder the Poynting 
flux is dissipated in the current sheets. The dissipation of Poynting flux is strongest at intermediate conductivities, 
a/rt ^1. At higher conductivities, the currents are more ideal and dissipationless, and at lower conductivities, 
dissipation is reduced because the conduction currents are weaker. Volume dissipation introduced in our resistive 
simulations implies that the light cylinder measurements of Poynting flux in Fig. |2] underestimate the true spin-down 
power of the resistive solutions, and should be treated as being within the error bands similar to the gray band of the 
force-free solution. 



